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Abstract Irrigation is the single largest anthropogenic 
water use, a modification of the land surface that signifi- 
cantly affects surface energy budgets, the water cycle, and 
climate. Irrigation, however, is typically not included in 
standard historical general circulation model (GCM) sim- 
ulations along with other anthropogenic and natural forc- 
ings. To investigate the importance of irrigation as an 
anthropogenic climate forcing, we conduct two 5-member 
ensemble GCM experiments. Both are setup identical to the 
historical forced (anthropogenic plus natural) scenario used 
in version 5 of the Coupled Model Intercomparison Project, 
but in one experiment we also add water to the land surface 
using a dataset of historically estimated irrigation rates. 
Irrigation has a negligible effect on the global average 
radiative balance at the top of the atmosphere, but causes 
significant cooling of global average surface air tempera- 
tures over land and dampens regional warming trends. This 
cooling is regionally focused and is especially strong in 
Western North America, the Mediterranean, the Middle 
East, and Asia. Irrigation enhances cloud cover and pre- 
cipitation in these same regions, except for summer in parts 
of Monsoon Asia, where irrigation causes a reduction in 
monsoon season precipitation. Irrigation cools the surface, 
reducing upward fluxes of longwave radiation (increasing 
net longwave), and increases cloud cover, enhancing 
shortwave reflection (reducing net shortwave). The relative 
magnitude of these two processes causes regional increases 
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(northern India) or decreases (Central Asia, China) in 
energy availability at the surface and top of the atmo- 
sphere. Despite these changes in net radiation, however, 
climate responses are due primarily to larger magnitude 
shifts in the Bowen ratio from sensible to latent heating. 
Irrigation impacts on temperature, precipitation, and other 
climate variables are regionally significant, even while 
other anthropogenic forcings (anthropogenic aerosols, 
greenhouse gases, etc.) dominate the long term climate 
evolution in the simulations. To better constrain the mag- 
nitude and uncertainties of irrigation-forced climate 
anomalies, irrigation should therefore be considered as 
another important anthropogenic climate forcing in the 
next generation of historical climate simulations and multi- 
model assessments. 

Keywords Climate modeling • Irrigation • Forcing • 
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1 Introduction 

General circulation model (GCM) simulations of the his- 
torical period (1850-present) use observation based esti- 
mates of natural (e.g., solar variability, volcanic eruptions) 
and anthropogenic [e.g., anthropogenic aerosols, green- 
house gas (GHG) concentrations, land use, etc.] forcings to 
simulate the last 150 years of climate evolution. Historical 
simulations provide an important baseline for understand- 
ing the context of future climate change, are used to test 
climate system responses to specific forcings, and are 
required for detection and attribution studies. These 
experiments therefore form an integral part of the model 
intercomparisons organized in tandem with the Intergov- 
ernmental Panel on Climate Change (IPCC) assessment 
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reports (e.g., Taylor et al. 2012). Recent historical model- 
ing studies have addressed diverse topics, including the 
impact of anthropogenic aerosol emissions on North 
Atlantic climate variability (Booth et al. 2012), the role of 
the ocean in the recent surface warming hiatus (Kosaka and 
Xie 2013), and the detection and attribution of anthropo- 
genically forced warming trends (Santer et al. 2013). One 
anthropogenic forcing that is typically not included in 
standard historical GCM simulations, including the IPCC 
model intercomparisons, is irrigation. 

Irrigation is the single largest anthropogenic water use, 
accounting for 70 % of global freshwater withdrawals and 
90 % of consumption (Siebert et al. 2010). Irrigation has 
expanded and intensified rapidly over the twentieth cen- 
tury, accelerating after 1950 with the increased exploitation 
of groundwater resources (Rodell et al. 2009; Scanlon et al. 
2012; Wada et al. 2010), and the climate effects of irri- 
gation have been broadly documented (e.g., in the review 
by Mahmood et al. 2013b). In areas where evapo transpi- 
ration rates are moisture limited, the addition of irrigation 
water shifts the surface energy balance to favor latent over 
sensible heating, decreasing the Bowen ratio and reducing 
soil and air temperatures (Kueppers et al. 2007 ; Sacks et al. 
2009). This first order effect has been documented pri- 
marily in areas characterized by growing seasons with low 
moisture availability and high atmospheric demand, such 
as the Central Valley of California (Bonfils and Lobell 
2007) and the southern Great Plains (Barnston and 
Schickedanz 1984). Enhanced evapotranspiration from 
irrigation can also increase atmospheric water vapor, with 
potential consequences for the regional and global radiative 
balance. For example, Boucher et al. (2004) found that 
irrigation, when added to a GCM, acted as a net positive 
radiative forcing due to the enhanced trapping of longwave 
radiation from increased humidity, an effect that may be 
enough to increase local minimum temperatures (Mah- 
mood et al. 2013a). Alternatively, if increased evapo- 
transpiration from irrigation also enhances cloud 
formation, this may reduce the net solar forcing at the 
surface or top of the atmosphere (TOA) by increasing 
albedo and shortwave reflection. Finally, although large 
uncertainties remain, there is evidence that irrigation may 
even enhance precipitation locally and downwind of major 
irrigation regions (e.g., Barnston and Schickedanz 1984; 
DeAngelis et al. 2010; Lo and Famiglietti 2013). 

There is clear evidence that irrigation has important 
local, regional, and possibly global effects on climate. 
Despite this, irrigation is typically not included as an 
anthropogenic forcing in historical GCM simulations along 
with land cover change, aerosols, GHGs, and other forc- 
ings. Here, we conduct a new suite of historical GCM 
experiments with the Goddard Institute for Space Studies 
(GISS) atmosphere GCM, ModelE2 (Schmidt et al. 2014), 


coupled to the Russell ocean model (Russell et al. 1995) 
and using historically estimated irrigation rates from the 
dataset of Wisser et al. (2010). Ours is the first GCM study 
of irrigation to quantify the importance of irrigation as a 
climate forcing in the context of other natural and 
anthropogenic forcings during the climate evolution of the 
twentieth century, using the same framework as the forced 
historical simulations from the Coupled Model Intercom- 
parison Project version 5 (CMIP5, Taylor et al. 2012). 

2 Data and methods 

2.1 GISS ModelE2 

All experiments are conducted with the coupled GCM 
ModelE2-R (2° x 2.5° latitude/longitude resolution), 
which includes the latest version of the state-of-the-art 
GISS atmosphere general circulation model (ModelE2), 
with updated physics for inclusion in the CMIP5 (Schmidt 
et al. 2014), and the Russell (R) ocean model (Hansen et al. 
2007; Russell et al. 1995). Historical simulations using 
current (Miller et al. 2014) and previous (Hansen et al. 
2007) versions of the model demonstrate the ability of 
ModelE2 to realistically reproduce historical responses to 
external climate forcings, including volcanic eruptions and 
anthropogenic GHGs. For our experiments, we use the 
version of model physics with non-interactive atmospheric 
composition (referred to as NINT in Miller et al. 2014), 
where concentrations of ozone and radiatively important 
aerosols are prescribed according to observation based 
decadal averages. 

The irrigation module in ModelE2-R uses historically 
varying, monthly irrigation rates from a 0.5° resolution 
reconstruction of twentieth century global hydrography 
(Wisser et al. 2010). Wisser et al. (2010) derive historical 
irrigation rates from 1901 to 2002 using the University of 
Frankfurt/FAO Global Maps of Irrigated Areas (Siebert 
et al. 2005a, b) (extended back to 1901 by Frey dank and 
Siebert 2008) and an offline model of terrestrial water 
balance and transport (WBMplus; Federer et al. 2003; 
Vorosmarty et al. 1998). WBMplus includes explicit rep- 
resentations of human activities that affect the water cycle: 
irrigation (Wisser et al. 2008) and reservoir operation 
(Wisser et al. 2010). For areas equipped for irrigation in 
WBMplus, an irrigation water flux is applied to the soil to 
compensate for climate observation forced soil moisture 
deficits, including an explicit treatment for paddy rice 
irrigation. In ModelE2-R, this new irrigation flux of water, 
calculated offline in WBMplus, is added to the top of the 
vegetated fraction of the soil column, below the vegetation 
canopy. Without a separate soil column for the irrigated 
crop fraction, irrigation is applied to the total vegetated 
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Fig. 1 Annual average 
irrigation rates (mm day -1 , log 
scale) for a 1901-1920 and 


b 1981-2000 



fraction. The implications of this simplification are 
addressed in the discussion section of this paper. Previous 
studies with the GISS model have used the irrigation 
module to investigate irrigation impacts on climate during 
the historical period (Puma and Cook 2010; Shukla et al. 
2013) and under increased GHG forcing (Cook et al. 201 1). 
Because our model simulations for this study begin in 
1850, we linearly extrapolated the model irrigation from 
1901 back to 1850. Further details on the development of 
the underlying irrigation dataset, and its implementation in 
ModelE2-R, can be found in Wisser et al. (2010) and Puma 
and Cook (2010), respectively. 

Irrigation is concentrated in the major regions of agri- 
cultural production in the Northern Hemisphere (Fig. 1). In 
the early part of the century (1901-1920, Fig. la), most 
irrigation occurs in Monsoon Asia, with secondary centers 
in the Mediterranean, the Middle East, and North America. 
By the end of the twentieth century (1981-2000, Fig. lb), 
irrigation rates intensify in these areas, and expand into 
additional regions, such as Northern Europe and Central 
Asia. Differences between the early and late twentieth 
century irrigation maps reflect the steady increase in global 
irrigation rates over the last 100 years (Fig. 2a), with a 
notable acceleration after about 1950 that coincides with 
the start of the green revolution and exploitation of 


groundwater resources. The vast majority of global irri- 
gation is concentrated in Monsoon Asia (Fig. 2b), where 
water intensive rice cultivation dominates the agricultural 
landscape. In North America there is a notable jump at 
1950, coinciding with the start of exploitation of the 
Ogallala aquifer (Scanlon et al. 2012). Over the European 
domain (Fig. 2d), irrigation increases over the course of the 
twentieth century before plateauing in the 1980s. 

2.2 Experiments 

We conduct two 5-member ensemble experiments with 
ModelE2-R from 1850-2000, with initial conditions taken 
from a near-equilibrium pre-industrial control run with 
atmospheric composition and solar irradiance held fixed at 
1850 values. In experiment HIST, all forcings are identical 
to those used in the standard ModelE2-R historical simu- 
lations submitted for the CMIP5 effort (Miller et al. 2014), 
including solar and volcanic forcing, anthropogenic GHGs, 
land cover, and tropospheric aerosols (natural and anthro- 
pogenic). In experiment HIST + IRR we include, as an 
additional forcing, transient historical irrigation rates, as 
described previously. Differences between the two exper- 
iments (HIST + IRR — HIST) are used to assess the cli- 
matic impact of irrigation. Because we are interested in the 


<£) Springer 


B. I. Cook et al. 


Annual Irrigation, km 3 




Fig. 2 Annual irrigation (km 3 ) from 1901-2000: a global, b Monsoon Asia (60°E— 130°E, 5°N— 40°N), c Western North America 
(125°W— 90° W, 30°N— 45°N), and d Europe and the Mediterranean (13°W— 40°E, 30°N— 70°N) 


forced response, all analyses are based on the ensemble 
mean from each model experiment. Statistical differences 
in the mean and median between the two simulations are 
assessed using the two-sided Student’s t test and Wilcoxon 
rank sum test, respectively. We focus our analyses pri- 
marily on the latter part of the twentieth century when 
irrigation rates, and their accompanying climate impacts, 
are largest. Unless otherwise noted, all anomalies (time 
series plots) are relative to the pre-industrial control run 
from which the initial conditions were taken and A refers to 
differences between the two experiments (HIST + IRR 
-HIST) for 1981-2000. 

3 Results 

3.1 Global average response 

In both experiments, the historical evolution of net 
radiation (Rnet) at the TOA follows the expected 
increase from enhanced GHG forcing and trapping of 


outgoing longwave radiation (Fig. 3a). Large magnitude, 
but short-lived, negative spikes in Rnet correspond to 
volcanic eruptions and the associated negative aerosol 
forcing. At the TOA, Rnet is about +0.02 W m : higher 
in HIST + IRR compared to HIST, similar to the mag- 
nitude of irrigation forcing (+0.03 W m 2 ) found in the 
experiments of Boucher et al. (2004). Irrigation forcing 
at the TOA in our experiments is insignificant 
( p > 0.05), however, suggesting that irrigation has a 
negligible effect on the global TOA radiative balance. It 
should be noted, though, that the standard radiative 
forcing concept might be difficult to apply to irrigation 
because of the temperature, energy balance, and water 
vapor feedbacks inherent in the irrigation-climate 
response (Boucher et al. 2004). 

Irrigation has a cooling effect on global average 
land + ocean surface air temperatures (Fig. 3b), but dif- 
ferences between the two simulations are small and insig- 
nificant {jp 0.05, AT mean — 0.07 K, AT me£ n an — 

—0.08 K). If the temperature average is calculated for land 
areas only, however, the cooling effect is larger and 
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Global Average Anomalies, Land+Ocean 




Fig. 3 Annual, globally (land and ocean) averaged anomalies of a net 
radiation at the top of the atmosphere (W m -2 ), b 2-m air 
temperature (K), c cloud cover (%), and d precipitation 
(mm day -1 ) for the ensemble averages of HIST and HIST + IRR. 

significant (A T mea „ = -0.17K ,p < 0.02;Zl7’ merfM „ = -0.l4K, 
p < 0.04) . This is a larger cooling effect than found in 
previous ‘realistic’ irrigation experiments. For example, 
land based cooling in Boucher et al. (2004) was only 
— 0.05K and previous irrigation experiments with ModelE 
using prescribed SSTs only caused cooling of — 0.10K 
(Puma and Cook 2010). Sacks et al. (2009) actually found a 
slight net annual warming (+0.015K), as dynamically 
induced warming during the winter balanced irrigation 
induced cooling during the summer (— 0.06K) and fall 
(— 0.04K). Cloud cover increases are highly significant in 
HIST + IRR (p< 0.001) over global land and ocean areas, 
with greater increases over land. Precipitation in HIS- 
T + IRR also increases over land and ocean in the irriga- 
tion run, but this increase is only marginally significant 

(p<0.10). 


Anomalies are relative to a 41 -year average from the pre-industrial 
control run used to generate the initial conditions for all ensemble 
members 


3.2 Temperature trends 

Increasing irrigation rates over the twentieth century cause 
significant differences in regional temperature trends. Both 
HIST (Fig. 4a) and HIST + IRR (Fig. 4b) show wide- 
spread and significant warming in annual average surface 
air temperatures from 1971 to 2000 (insignificant differ- 
ences at p > 0.05 in Fig. 4 have been masked out). In 
many areas with major irrigation expansion or intensifica- 
tion, however, these warming trends weaken (Western 
North America), become insignificant (Europe), or, in 
some cases, reverse to cooling trends (extra-tropical South 
America, India). 

The time when irrigation begins to significantly affect 
surface air temperature varies by region. In Central Asia 
and China (Fig. 5a, b), temperatures in the two simulations 
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Fig. 4 Ensemble average 
annual surface air temperature 
trends (1971-2000, K year' 1 ) 
from the a HIST and b 
HIST + IRR simulations. 
Insignificant trends (p > 0.05) 
have been masked out 



begin to diverge after about 1950. Both regions show 
warming up to the 1950s, followed by multi-decadal 
cooling often ascribed to increased aerosol forcing (e.g.. He 
et al. 2013; Qian and Giorgi 2000). Despite the temperature 
offset, both regions have similar temperature trends over 
1971-2000. This contrasts sharply with India (Fig. 5c), 
where irrigation rates are high enough across the whole 
century to lead to continuously cooler surface temperatures 
in HIST + IRR. Temperatures trends do diverge over 
India, with increasing irrigation causing a slight cooling 
trend over the last 30 years of the century. Late twentieth 
century temperature trends are similar between the HIST 
and HIST + IRR simulations for the Mediterranean 
(Fig. 5d), the Middle East (Fig. 5e), and Western North 
America (Fig. 5f). Cooler average temperatures in HIS- 
T + IRR manifest in the Mediterranean and Western North 
America beginning in the 1970s and in the 1930s in the 
Middle East. 

3.3 Spatial response: 1981-2000 

Irrigation effects on climate in the late twentieth century 
(1981-2000) are highly localized, coinciding with land 


areas where irrigation is most intense. Irrigation causes 
widespread cooling (Fig. 6a), primarily over Western 
North America, the Mediterranean, the Middle East, India, 
Indochina, Central Asia, and China. Localized cooling is 
generally on the order of —0.5 to — 1.0K, with much more 
intensive cooling over northern India, where the highest 
irrigation rates in the model and irrigation dataset occur. 
Irrigation increases atmospheric moisture, primarily at low 
levels, and this is reflected in the 850 hPa specific humidity 
(Fig. 6b) and low level (surface to 680 hPa) cloud cover 
(Fig. 6c). In concert with these humidity increases, pre- 
cipitation is also enhanced (Fig. 6d), especially in Western 
North America, the Middle East, and Central Asia. The 
only major region where annual land precipitation signifi- 
cantly declines is in far eastern India, Bangladesh, and 
Myanmar in Monsoon Asia due to drying during the 
summer monsoon. 

Despite little impact on the globally averaged radiative 
balance, irrigation has a large effect on regional radiative 
fluxes (Fig. 7). Significant changes in net total radiation at 
the surface (Fig. 7a) and TOA (Fig. 7b) occur primarily 
over Monsoon and Central Asia. Over the Indian subcon- 
tinent and Southeast Asia, net radiation increases, while 
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2-m Air Temperature Anomalies (K) 



2000 



Fig. 5 Regional ensemble average annual surface temperature 
anomalies (K) from major irrigation regions: a Central Asia 
(72°E~98°E, 35°N-45°N), b China (100°E-130°E, 20°N-45°N), 
c India (70°E— 95°E, 5°N— 35°N), d the Mediterranean 

(13°W— 30°E, 30°N— 45°N), e Middle East (31°E-68°E, 10°N-44° 
N), and f Western North America (125°W-90°W, 30°N-45°N). 


Temperature curves are smoothed with a 5-year lowess filter to 
emphasize long term trends and variability. Black dashed lines are the 
best fit linear regression lines from 1971 to 2000, based on a fit to the 
unsmoothed ensemble average temperatures. Anomalies are relative 
to a 41 -year average from the pre-industrial control run used to 
generate the initial conditions for all ensemble members 


over Central Asia and China net radiation decreases. Irri- 
gation reduces net shortwave at TOA (Fig. 7c) and the 
surface (Fig. 7d), especially over Western North America, 
the Middle East, and Asia. These reductions coincide clo- 
sely with the increases in cloud cover (Fig. 6c), indicating 
this change is through increased reflection by clouds. 
Indeed, areas where net shortwave is actually increasing 
(e.g., northern India) are areas where cloud cover declines. 
Irrigation increases net longwave fluxes (Fig. 7e, f). In 
ModelE2-R this occurs because the cooler soil tempera- 
tures reduce outgoing longwave fluxes; increased longwave 
trapping from enhanced near surface humidity is small (not 
shown). Over North America, the longwave and shortwave 
effects largely cancel, leading to negligible changes in net 
radiation. In India, the enhanced longwave dominates, 
increasing net total radiation, while in Central Asia and 
China, the shortwave reductions drive an overall decline in 
energy availability. 


Irrigation has a much larger impact on surface energy 
partitioning (sensible, Fig. 8a, and latent. Fig. 8b) than on 
total energy availability in our experiments. In areas where 
irrigation rates are high (Western North America, the 
Mediterranean, Middle East, India), the addition of irriga- 
tion water shifts the Bowen ratio to favor latent over sen- 
sible heating, and these shifts coincide closely with much 
of the surface cooling in Fig. 6a. Outside the main irriga- 
tion regions, significant increases in latent heating and 
reductions in sensible heating are related to local increases 
in precipitation (e.g., the Sahara; Fig. 6d). The dominance 
of Bowen ratio versus energy balance shifts for the tem- 
perature response in our experiments contrasts sharply with 
the irrigation study of Sacks et al. (2009). In their study. 
Sacks et al found instead that the cooling response was 
driven primarily by irrigation-induced increases in cloud 
cover and reductions in surface shortwave radiation, rather 
than by direct evaporative cooling. 
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ALow Cloud Cover (%) 



APrecipitation (mm day 1 ) 



Fig. 6 Irrigation effects (HIST + IRR — HIST) on a 2-m surface air 
temperature, b specific humidity at 850 hPa, c low cloud cover 
(surface to 680 hPa), and d precipitation. Differences are based in the 


ensemble means of each experiment for 1981-2000. Insignificant 
differences (p > 0.05, two-sided Student's t test) have been masked 
out 


3.4 Regional and seasonal effects 

Highest irrigation rates occur in Asia, and this is where the 
largest temperature and precipitation responses occur 
(Fig. 9). In the relatively temperate mid-latitude regions of 
China (Fig. 9, top row) and Central Asia (Fig. 9, middle 
row), irrigation is highest during boreal summer, and it is 
during the warm months that most of the significant cool- 
ing occurs (indicated by red stars, p < 0.05, two-sided 
Student's t test). Irrigation significantly enhances warm 
season precipitation in June in China and during May, 
June, August, and September in Central Asia. Irrigation 
rates are much higher over the region of India + Bangla- 
desh, with larger (and significant) magnitude temperature 
and precipitation anomalies in most months (Fig. 9, bottom 
row). This is partly due to the more sub-tropical nature of 
this region, characterized by relatively high evaporative 
demand throughout the year and a pronounced dry season 
during boreal winter. Most cooling in India + Bangladesh 
occurs during this dry season, but precipitation responses 
are seasonally dependent. Precipitation increases are small, 
but significant, during the dry season, and are likely linked 
to increased evapotranspiration from irrigation during the 
relatively warm dry season. During the main summer 


monsoon season, precipitation deficits are large and con- 
centrated in eastern India and Bangladesh. These changes 
may be related to the tendency for irrigation to weaken the 
summer monsoon circulation in this region (discussed in 
detail in Shukla et al. 2013). An analysis of the surface 
energy balance partitioning in these regions (Fig. 10) 
indicates that the temperature changes are due primarily to 
changes in the Bowen ratio, rather than total energy 
availability. In China and India + Bangladesh, significant 
changes in net radiation at the surface (Fig. 10, left col- 
umn) are all positive, indicating the dominance of the 
longwave effects, while in Central Asia increased short- 
wave reflection causes a net reduction in surface energy. 
These changes in net energy availability, however, are 
much smaller than the changes in latent (Fig. 10, center 
column) and sensible (Fig. 10, right column) heat fluxes. 

Irrigation also has a large and significant impact on 
climate in the Mediterranean (Fig. 1 1 , top row). Middle 
East (Fig. 11, middle row), and Western North America 
(Fig. 1 1 , bottom row). These three regions are generally 
characterized by warm summers with high evaporative 
demand and, as with China and Central Asia, the warm 
season is when the temperature (cooling) and precipitation 
(enhanced) responses are largest. Notably, precipitation 
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Fig. 7 Same as Fig. 6, but for net total radiation at the a top of the atmosphere and b surface, net shortwave radiation at the c top of the 
atmosphere and d surface, and net longwave radiation at the e top of the atmosphere and f surface 


increases in Western North America occur during July- 
August-September, the same season of enhancement found 
in the irrigation experiments of Lo and Famiglietti (2013). 
Irrigation causes some net increases in net radiation 
(Fig. 12, left column) in the Mediterranean and Middle 
East, but the shortwave and longwave effects cancel out 
over Western North America. As with Asia, shifts in latent 
(Fig. 12, center column) and sensible (Fig. 12, right col- 
umn) heat fluxes are much larger than any changes in net 
energy availability. 

4 Discussion and conclusions 

Human land use and land cover change is an important 
anthropogenic climate forcing (Mahmood et al. 2013b), 
with regional impacts that may match or even exceed 
localized responses to globally well-mixed anthropogenic 


GHG forcing (Pielke et al. 2011). But while the effects of 
land cover change on historical and future climate have 
been comprehensively assessed (e.g., Feddema et al. 2005; 
Kumar et al. 2013; Pitman et al. 2009), only recently have 
studies begun to analyze the response of the climate system 
to second order land use effects (Luyssaert et al. 2014). 
This includes cropland management (Lobell et al. 2006), 
crop growth and development (Chen and Xie 2012), and 
irrigation (Lobell et al. 2006; Sacks et al. 2009). In this 
study, we use a new set of historically forced GCM 
experiments to investigate the importance of irrigation as 
an historical climate forcing. 

Compared to other regionally important anthropogenic 
forcings, such as land cover change (Boisier et al. 2012) 
or aerosols (e.g., Bellouin et al. 2013; Booth et al. 2012), 
irrigation has relatively small effects on the total energy 
availability at the surface or TOA for most areas. Sig- 
nificant perturbations to net total radiation are largely 
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Fig. 8 Same as Fig. 6, but for 
a surface sensible heat flux and 
b surface latent heat flux 
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confined to Asia (Fig. 7a,b), and range between about —4 
(e.g., Central Asia and China at TOA) to +10 (e.g., India 
at the surface) W m 2 , depending on whether shortwave 
or longwave effects dominate. Instead, irrigation has a 
much larger and widespread influence on the surface 
partitioning between sensible and latent heat fluxes 
(Fig. 8). It is these shifts in the Bowen ratio that cause 
widespread regional cooling and enhancement of pre- 
cipitation, with climate responses generally scaling with 
the intensity of irrigation. The small impact of irrigation 
on net energy availability is perhaps not surprising, given 
that the addition of irrigation does not directly add 
substantial amounts of energy to the global earth system. 
Any changes in net radiation must then come from 
secondary effects (e.g., cloud cover changes, increased 


atmospheric water vapor, etc.) which, as noted, were 
small at the global scale. 

Contrary to observations (e.g., Jones et al. 2012), 
warming trends in HIST + IRR from 1971 to 2000 become 
insignificant or actually reverse in Northern Europe and 
India, suggesting that climate in ModelE2-R is overly 
sensitive to irrigation, at least in certain regions. This 
overestimation of irrigation effects on climate, and the 
subsequent divergence between our GCM simulations and 
observed warming trends, could occur for a variety of 
reasons. First, the magnitude of climate responses to irri- 
gation depends crucially on the evaporative regime, i.e. the 
degree to which evapotranspiration rates are energy 
(demand) versus moisture (supply) limited (Cook et al. 
2011; Puma and Cook 2010). The extent to which models 
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Fig. 9 For China (100°E-130°E, 20°N-45°N, top row). Central 
Asia (72°E--98°E, 35°N— 45°N, middle row), and India 
(70°E— 95°E, 5°N— 35°N, bottom row), monthly irrigation rates in 
the HIST + IRR ensemble (left column), and monthly differences 
(HIST + IRR — HIST) in 2-m air temperature ( center column) and 


precipitation (right column). Differences are calculated using the 
ensemble averages for 1981-2000. Significant differences at the 
p < 0.05 level, based on a two-sided Student’s t test, are indicated by 
the red stars 


can accurately simulate observed evaporative regimes, 
however, is affected by numerous model dependent factors, 
including precipitation biases and land surface parameter- 
izations (e.g., rooting depth, vegetation seasonality, etc.). 
Ultimately, if model soils are too dry compared to the real 
world, irrigation will have a larger effect on evapotrans- 
piration and the subsequent climate response will be 
amplified. Relatedly, land-atmosphere coupling strengths 
are known to vary widely across GCMs (Koster et al. 
2006). Second, models with lower resolution than the 


original irrigation dataset (0.5°) may amplify irrigation- 
climate responses (Puma and Cook 2010; Sacks et al. 
2009). For this study, the underlying irrigation dataset was 
aggregated to the coarser resolution (2° x 2.5°) of Mo- 
delE2-R by summing all irrigation grid cells that would fall 
within each coarser ModelE2-R grid cell. A single cell at 
coarser resolution may then include both irrigated and non- 
irrigated areas from the original irrigation dataset, even as 
the whole cell is treated as irrigated in the model. This can 
result in greater irrigated areas and potentially higher 
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Fig. 10 For the same regions in Fig. 9 (China, Central Asia, India), differences (HIST + IRR — HIST) in net radiation at the surface (left 
column), surface latent heat flux ( center column ), and surface sensible heat flux (right column ) 


evapotranspiration rates and climate responses. To improve 
the realism of these responses, increases in spatial resolu- 
tion of the model’s land and vegetation modules would 
reduce overestimates of irrigated area. More directly, 
ModelE2-R should be modified to simulate the dynamics 
of a separate soil column for irrigated crops within each 
grid cell. These changes would lead (ideally) to improved 
soil moisture and evapotranspiration dynamics for both 
irrigated crops and other vegetation. Finally, there are still 
large uncertainties in how irrigation interacts with other 
anthropogenic forcings in GCMs. For example, anthopo- 
genic aerosol forcing has important effects on the surface 
energy balance in many regions, which may affect the 


evaporative regime and evapotranspiration rates (Naza- 
renko and Menon 2005). 

Where possible, we attempted to compare our results 
with other irrigation studies, but a comprehensive com- 
parison was difficult for multiple reasons. For example, 
there is little standardization across modeling studies in 
the experimental setup, amount of irrigation forcing 
applied, and the model diagnostics reported. Analyses of 
irrigation effects on climate would therefore clearly 
benefit from standardized multi-model intercomparisons, 
similar to existing projects studying other land surface- 
climate interactions, such as GLACE (Koster et al. 2006; 
Seneviratne et al. 2013) and LUCID (Boisier et al. 2012; 
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Fig. 11 For the Mediterranean (13°W— 30°E, 30°N— 45°N, top row), 
Middle East (31°E— 68°E, 10°N— 44°N, middle row), and Western 
North America (125°W— 90°W, 30°N— 45°N, bottom row), monthly 
irrigation rates in the HIST + IRR ensemble (left column), and 
monthly differences (HIST + IRR — HIST) in 2-m air temperature 


(center column) and precipitation (right column). Differences are 
calculated using the ensemble averages for 1981-2000. Significant 
differences at the p < 0.05 level, based on a two-sided Student's t test, 
are indicated by the red stars 


Pitman et al. 2009). For irrigation, these intercompari- 
sons are becoming increasingly tractable as GCM groups 
have begun integrating their own irrigation schemes, 
including the Community Earth System Model (CESM, 
Levis and Sacks 2011) and Laboratoire de Meteorologie 
Dynamique (LMD-Z, Guimberteau et al. 2012). Ideally, 
such an intercomparison would, at minimum, (1) stan- 
dardize the set of model diagnostics reported and ana- 
lyzed (e.g., energy balance terms, temperature, etc.) and 
(2) use a common dataset of irrigation forcing, such as 
the WBMplus data set used here. Such a set of experi- 
ments would be especially helpful for highlighting model 


based uncertainties in the irrigation-climate response 
(e.g., model resolution, surface climate biases, etc.). 
Additionally, given the broad evidence presented here, 
and in other modeling and empirical studies, for a sig- 
nificant impact of irrigation on climate, we believe that 
irrigation should be included as another standard histor- 
ical anthropogenic climate forcing in future multi-model 
assessments, such as the CMIP6 (Meehl et al. 2014). The 
importance of irrigation for historical (and future) cli- 
mate can then be assessed in comparison to other 
important anthropogenic forcings, such as aerosols, land 
cover change, and GHGs. 
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